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THEORETICAL STUDY OF HULL-ROTOR AERODYNAMIC 
INTERFERENCE ON SEMIBUOYANT VEHICLES 

by 

S. B. Spangler and C. A. Smith 
Nielsen Engineering & Research, Inc. 

SUMMARY 

Analytical methods have been developed to predict the pressure 
distribution and overall loads on the hulls of airships which have 
close coupled, relatively large and/or high disk loading propul- 
sors for attitude control, station keeping, and partial support of 
total weight as well as provision of thrust in cruise. The methods 
comprise a surface-singularity, potential-flow model for the hull 
and lifting surfaces (such as tails) and a rotor model which cal- 
culates the velocity induced by the rotor and its wake at points 
adjacent to the wake. Use of these two models provides an invis- 
cid pressure distribution on the hull with rotor interference. A 
boundarylayer separation prediction method is used to locate 
separation on the hull, and a wake pressure is imposed on the 
separated region for purposes of calculating hull loads. 

A number of calculations were made to investigate and illus- 
trate the nature of hull-rotor interactions. One basic configura- 
tion was examined for hull alone, rotor alone, and hull-rotor 
combinations in hover, axial flow, and crossflow (8-90°) in and 
out of ground effect. Comparison with the limited data available 
on this configuration showed generally good agreement on the hull 
alone, including cases of large separated flow on the hull, and 
the rotor alone. The predicted loads on the hull-rotor combina- 
tions appear reasonable and consistent, but differ from measured 
values in sign of the normal force. The quality and quantity of 
the data are insufficient to establish the source of the differ- 


ences. 


INTRODUCTION 


Historically, airships have been rather lightly powered and 
have used their propulsion systems to overcome drag in cruise. 

The airships were operated close to the neutrally buoyant condi- 
tion, and aerodynamic lift on the hull was used to make up differ- 
ences between buoyant force and weight. 

More recently, interest has been generated in airships which 
have a considerably higher installed thrust capability relative 
to vehicle weight. Recent NASA-sponsored work has identified the 
Heavy Lift Airship (HLA) and a passenger feeder liner as promis- 
ing applications for modern airship technology. These vehicles 
characteristically have relatively large rotors or propellers, 
which can be mounted at various locations on the hull and which 
are expected to produce thrust for attitude control, station 
keeping, and partial support of the total weight as well as for 
overcoming drag in cruise. Because of the higher thrust levels 
produced and the tilting capability in these propulsors, there is 
the possibility of large interactions that can occur between the 
rotor wakes and the hull, which can produce either favorable or 
unfavorable interference loads on the vehicle. Due to the mag- 
nitudes of these effects, it is necessary to know their charac- 
teristics before the design of a vehicle is finalized and to use 
their characteristics to determine the control power available, 
the capability to station-keep in a wind, and the approach to the 
control system design. 

Since the problem is relatively new, there are no existing 
analytical methods to predict interference effects of this nature 
and very little appropriate data. Perhaps the most applicable 
data source is results of tests conducted under Goodyear Aerospace 
Corporation sponsorship at the Ames Research Center by Nielsen 
Engineering & Research on a l/79th scale model of the Heavy Lift 
Airship vehicle being studied by GAC under their NASA contract.^ 
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The objective of the program described in this report is to 
develop a theoretical method for predicting hull loads with rotor 
interference on semibuoyant airships. The method is to be general 
in the sense that it can be applied to a variety of hull shapes 
(not necessarily axisymmetric) with rotors, of a wide range of disk 
loading, located adjacent to the hull. In addition, it is neces- 
sary that the method be applicable to a vehicle operating near 
the ground and to vehicle attitudes relative to the wind direction 
such that substantial areas of separated flow may occur on the hull. 

This report describes the development of the methods and re- 
sults obtained with the method on typical configurations. In- 
cluded are comparisons with experimental data whenever possible. 
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NOMENCLATURE 


b 


Number of rotor blades 



C X' C Y' C Z 


D 


Moment coefficients on hull* nondimensionalized 
by q.V 

Pressure coefficient* based on 

Modified pressure coefficient* based on U * 
equation 16 

Rotor thrust coefficient* T/p^ (fiR) 2 J1R 2 

Force coefficients on hull* nondimensionalized 
by q.V 2 /'* 

hull maximum diameter 


F{s) 

H 

L 


“t 


M x' M y' M z 


n 


Empirical parameter in separation criterion 

Height of hull geometric center above ground 

Coordinate along axis of rotor wake 

Rotor thrust moment, or moment about the flapping 
hinge due to the thrust 

Moments on hull about x, y* and z axes* respec- 
tively, positive "right wing down", nose up, 
and nose right, respectively (fig. 1) 

Coordinate normal to hull surface 


r 

-+ 

r 


Free stream dynamic pressure 
Radial coordinate 

Vectorial distance from origin to a field point 



Re 


s 

s' 


m 


s 


Rotor radius 

Distance from point on rotor disk to field point 

Free stream Reynolds number, based on hull 
maximum diameter 

Reynolds number based on distance along hull 
surface from virtual origin 

Distance from rotor center to field point 

Coordinate along hull surface in a y,z plane 

Distance along hull surface from stagnation point 
to location of virtual origin in separation 
criteria 

Distance along hull surface from stagnation point 
to location of maximum velocity 

Vector from rotor center to point on wake 
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NOMENCLATURE (continued) 


Rotor disk area 


T 


Rotor thrust 


u,v,w 


u 


e 


U 


m 


Velocity components in x,y,z directions, 
respectively 

Magnitude of inviscid (or boundary layer edge) 
velocity 

Magnitude of maximum inviscid velocity on hull 


U„ 


Tangential component of relative velocity at 
rotor blade element 


U a'V W a 


U c' V c' W c 


u ,v ,w 

r r i 


U R' V R' W R 


V 

V 
V. 


Components of rotor- induced velocity due to 
axial vorticity 

Components of rotor-induced velocity due to 
circumferential vorticity 

Components of rotor-induced velocity due to 
radial vorticity 

Components of rotor-induced velocity due to 
root vorticity 

Hull volume 

Total velocity 

Rotor wake velocity at plane of rotor (from 
momentum theory) 

Rotor- induced velocity due to root vorticity 

Free stream velocity 


x,y,z 


X,Y,Z 


V y Q' Z Q 

W Z R 


a 


Coordinates in a body axis system centered at 
the body geometric center; x forward, y to 
starboard and z downward (fig. 1) 

Forces on hull along the x,y, and z directions, 
respectively 

Coordinates of rotor center in hull coordinate 
system (fig. 1) 

Coordinate system with origin at center of rotor 
hub; x R downstream, y R to starboard and z R 
vertically upward (fig. 1) 

Angle of attack 
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NOMENCLATURE (concluded) 


Angle of eldeellp 

Rotor blade circulation strengths due to, respec- 
tively, root vorticity, ring vorticity, axial 
vorticity, and radial vorticity (fig. 2) 

Total circulation on rotor blade 

Part of blade circulation independent of azimuth 

Portion of blade circulation that varies with the 
azimuthal angle 4* 

Circulation strength of the rotor blade axial 
vorticity 

Circulation strength of the rotor blade circum- 
ferential vorticity 

Circulation strength of the rotor blade root 
vorticity 

Strength of a vortex filament 

Rotor inflow ratio, nondimen sional component of 
velocity perpendicular to rotor dish, 

V„ sin o + V. 

_Z I 

ft R 

Rotor tip-speed ratio, nondimensional component 

of velocity parallel to rotor disk, 

V_ cos o 
00 

Hr 

Free- stream density 

Surface source density distribution 
Velocity potential for disturbance velocity 
Velocity potential for total velocity 
Rotor skew angle, figure 1 

Azimuthal angle around the rotor disk, figure 1 

Rotor rotational speed 

Magnitude of rotor axial vorticity 

MagnituJ ' of rotor circumferential vorticity 

Rotor tip vortex 


OVERALL APPROACH 


The overall approach to the development of a hall-rotor in- 
terference model is to use potential flow methods to describe the 
flow over a hull and the flow field induced by a rotor, combine 
these two methods to calculate an inviscid hull pressure distri- 
bution with hull-rotor interference, and modify this pressure 
distribution according to predicted flow separation and wake 
pressure calculations. Thus, two potential flow models, a 
method for predicting boundary layer separation, and a wake 
pressure model are required. 

The hull potential flow method used is the incompressible 
lifting model of reference 4. This method, and the resulting 
computer program (reference 5) , calculates the inviscid pressure 
distribution on an arbitrary nonlifting body, a lifting wing, or 
a body-wing combination. The program places the configuration 
in a uniform onset flow and uses various types of singularities, 
depending on the surface being modeled. It has a provision for 
symmetry planes to reduce computation time. For purposes of 
this investigation, the program can be made to model an airship 
hull of arbitrary shape, in or out of ground effect, with or 
without tail surfaces. 

The rotor potential flow model is based on the theory devel- 
oped by Heyson (refs. 6 and 7) for a rotor alone in a uniform 
flow. The model uses vortex filaments to represent the blade 
bound, trailing, and shed vorticity and calculates the strength 
of the vorticity distribution from the rotor thrust, advance ratio 
and angle of attack. The vorticity distribution can be used to 
calculate rotor-induced velocities in the vicinity of the rotor. 

When the two models are combined, the influence of the hull 
on rotor thrust can be estimated, the rotor vorticity distribu- 
tion calculated, the rotor wake placed in relation to the hull, 
and the rotor-induced velocities used as a nonuniform onset flow 
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added to the uniform fre*-stream flow to impinge on the hull. The 
hull-tail surface singularity distribution is then calculated in 
the presence of the two onset flows to obtain the inviscid pres- 
sure distribution. 

Boundary- layer separation location is computed by applying 
empirically developed separation criteria to the hull inviscid 
pressure distribution. For those parts of the hull that are lo- 
cated downstream of separation, the inviscid pressure is replaced 
by a wake pressure assumed constant over the separated region at 
the value at separation. Loads on the hull are then computed. 

The individual methods and their use in combination are des- 
cribed in following sections. 

FLOW MODELS 

Hull Potential Flow Model 

General Theory .- The hull flow model determines the forces 
and moments acting on the hull and the pressure distribution on 
the hull surface by representing the hull surface by a finite 
number of surface source singularities. The flow field used to 
calculate the required strengths of these source singularities in- 
cludes both the uniform free stream and the nonuniform flow in- 
duced on the hull surface by the rotor wakes. 

The mathematical solution to this problem proceeds in the 
following manner. At any point the fluid velocity is assumed to 
satisfy Laplace's equation which can be expressed as 

V 2 * = 0 (1) 

where 

V = grad * (2) 
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In addition the velocity is subject to two boundary conditions. 

The first is that the normal velocity component on the surface of 
body is zero and the other is that the velocity approaches that 
of the free stream at infinity. Symbolically, these boundary con- 
ditions are written 

ii « 0 (3) 

3n surface 


lim (grad $) - (4) 

> 

r+“ 

Divide the velocity potential into a part due to the uniform free 
stream and a part due to the disturbance velocity such that 

* 4*^ + 4* (5) 

where 

grad 4\ * V ro (6) 


The governing equation and boundary conditions then become 

V 1 * - 0 (7) 


3 ft 

3n surface 


4 + 

n • V c 

surface 


( 8 ) 


lim (grad 4) » 0 (9) 

*■ 

j* *»\> 

Use is now made of the result shown in reference 8 that the body 
surface can be covered by a surface source density distribution 
o. In terms of this source distribution the governing equation 
and boundary conditions given by equations (7) , (8) , and (9) re- 

duce to the single equation 
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( 10 ) 


2ito (p) 

surface 

where points p and q are on the surface of the body and dS 
is an elemental surface area of the body. This equation is a two- 
dimensional Fredholm integral equation of the second kind over 
the surface. 

Solution .- There are several methods to solve equation (10) , 
The one chosen for the present application is that developed at 
Douglas by Hess and Smith and reported in references 4 and 9 . 

This method consists of approximating the body surface by a large 
number of plane quadrilaterals. The source density is taken to 
be constant over each quadrilateral or panel and one point is 
specified (control point) on each panel where equation (10) is re- 
quired to hold. This reduces the integral equation to a set of 
linear equations equal to the number of panels used to represent 
the body surface. 

The resulting set of linear equations is used to determine 
the source densities. Once this is determined the flow veloci- 
ties at the control points are evaluated by calculating the con- 
tribution of each of the source elements and adding the components 
of the free stream. Finally, the pressure distribution over the 
surface is determined with the Bernoulli pressure equation, i.e. 

C p - 1 - HI) 

* OO 


o (q) dS ■ -n(p) *V 


- r 


Once the pressure distribution over the surface is determined the 
forces and moments are evaluated by assuming the pressure at a 
panel to act on the surface area of that panel to provide a force 
along the normal to the panel at the control point and summing 
over the body surface. The details of the computation and pro- 
gram are presented in reference 5. 
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The potential flow computer program developed at Douglas was 
modified in two ways in order to study rotor-hull configurations. 
The first addition was to compute the velocity at each control 
point induced by the rotor wakes and then add this velocity to 
the free-stream velocity to get the total onset flow. The flow 
model used to evaluate these rotor-wake-induced velocities is 
described in the following section. The second modification takes 
the inviscid pressure distribution determined by equation (11) 
and uses it to calculate the point of turbulent separation. 

The details of the flow separation model are described in 
a subsequent section. All panels downstream of the separation 
point are then assigned a pressure equal to the pressure at sep- 
aration. This modified pressure distribution is then used to 
calculate forces and moments. 

Lifting Surface Model 

The Douglas program is capable of representing body-lifting 
surface combinations in addition to bodies alone. In the present 
application, this capability provides a means of modeling tail 
surfaces on airship hulls, for example. The general approach is 
similar to that for the body alone. 

The lifting surface is divided into an arrangement of panels 
on the pressure and suction surfaces. A surface source distri- 
bution on both surfaces is used to model the thickness. The lift 
is modeled by a surface vorticity distribution, which is treated 
as an equivalent dipole distribution for calculation purposes, 
and trailing vorticity. The chordwise strength distribution of 
the surface vorticity at a given span station is uniform, and its 
strength is calculated from a Kutta condition. Reference 4 gives 
considerable information on the method. Any lifting surface plan- 
form, thickness distribution, and camber can be modeled by the 
method. The lifting surface may be placed anywhere on the body 
surface. The junction between the body and lifting surface is 
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handled by extending the surface one row of panels inside the 
body to eliminate the discontinuity in the lifting surface vor- 
ticity representation at the junction. 

Rotor Wake Model 

A potential flow model of a rotor wake is required to cal- 
culate the disturbance velocity induced by that rotor at points 
adjacent to and outside of the wake (i.e., on the hull). Gener- 
ally, the propulsors that are being considered for semibuoyant 
airships are low solidity devices like propellers and helicopter 
rotors rather than high solidity devices like turbofans. Conse- 
quently a rotor model rather than a jet model seems most appro- 
priate for the propulsor wake. 

The model chosen to represent the rotor wake is an extension 
of the model developed by Heyson (refs. 6 and 7). This model has 
the advantage of being an analytical model which represents most 
of the important features of the wake vorticity distribution and 
which is very quick to execute on a computer. The alternative is 
one of several numerical models which have been developed more 
recently and which can calculate in greater detail the rotor disk 
inflow distribution and velocity distribution within the wake, 
but at a great increase in computer time. Since the present 
interest is in predicting velocities outside the wake, the Heyson 
model is considered the best choice. 

The Heyson model is shown in figure 2. The circulation 
along the blade is taken to be uniform in the radial direction. 
Thus, trailing vortices emanate from the blade tip and the blade 
root. The root vortex (Yj) is a single vortex filament and lies 
along the wake centerline while the tip vortex forms the wake 
outer layer and is in the shape of a helix. The wake extends 
downstream to infinity along a straight centerline and the wake 
cross section is constant in size. 
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It is now assumed that the tip vortices are so closely spaced 
that they form a continuous cylindrical sheet of vorticity. This 
assumption is equivalent to assuming an infinite number of blades 
and constrains the analysis to determining the time-averaged in- 
duced velocities. The cylindrical vortex sheet representing the 
tip vortex is then divided into two components: one part being 

circumferential vorticity (or vortex rings) (y 2 ) and the other part 
being axial vorticity (y 3 ) * Hence the tip vortex is separated 
into two parts, the sum of which is the strength of the tip vor- 
tex. This is shown in figure 3 taken from reference 7. 

A rotating blade in forward flight will experience a change 
in the velocity field with a change in azimuthal position around 
the rotor disk. To account for this the blade circulation (and, 
therefore, the wake vorticity) is allowed to vary with the azimuth 
angle, ip, as sin ip. Thus additional radial vorticity (y 4 ) is 
shed from the trailing edge of the blades as a result of the 
changes in blade circulation. This change in circulation is con- 
tinuous, so the wake cylinder is completely filled with radial 
vorticity of a strength proportional to the derivative of the 
circulation change with azimuth. The total blade circulation is, 
therefore, given by the relation 

r = r Q + r sin n> (12) 

Further assumptions of the model are that the wake is non- 
contracting, is straight, and leaves the rotor disk in the pres- 
ence of a free stream at an angle (skew angle) determined by 
momentum theory. Finally, it is assumed that the wake does not 
impinge on any nearby surface nor is it influenced by either near- 
by surfaces or a non-uniform free stream. 

When determining the contribution to the induced velocity 
from each of the various components of vorticity, no simplifica- 
tions have been made as in the original development by Heyson 
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(ref. 7) where the axial vorticity was assumed negligible for all 
but high speed forward flight. 


The velocity at a point induced by a vortex filament is 
given by the Biot-Savart law. For a single filament this is ex- 
pressed as 


V (r) 



dtx(r-s) 



3 


(13) 


where F v is the strength of the filament. In the present situa- 
tion the induced velocity at a point is determined by the inde- 
pendent application of the Biot-Savart law to each of the four 
vortex systems in the wake. The total induced velocity is then 
the summation of these four contributions. 

The expression for the velocity due to the root vortex can 
be solved directly in closed form. However, the expressions for 
the other three velocity contributions must be solved by numeri- 
cally integrating around the rotor disk. These expressions are 
derived and presented in Appendix A. 

The application of the method requires specification of the 
thrust coefficient, the rotor radius, the angular rate, the free 
stream velocity and its angles of attack and sideslip, and the 
attitude of the rotor relative to the hull. From this informa- 
tion, the rotor thrust and the blade circulation can be calcula- 
ted assuming that the thrust produces a uniform downwash through 
the rotor disk in accordance with the momentum equations. Final- 
ly, the wake skew angle is calculated as the resultant direction 
of the vector sum of the rotor- induced downwash velocity and the 
horizontal and vertical components of the free-stream velocity. 
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HULL FLOW SEPARATION 


Separation Location 

The flow over the hull at all practical wind speeds will 
separate and change the hull pressure distribution from that pre- 
dicted by potential theory. In a largely axial flow condition, 
separation will occur over the aft end of the hull and the pres- 
ence of the rotors will have little influence on hull flow sepa- 
ration (unless there is a stern-mounted rotor) , nor will hull 
flow separation have much influence on the rotors. In a largely 
crossflow condition, separation occurs along the length of the 
hull, and the hull flow separation - rotor performance interaction 
can be significant. A flow separation model was incorporated in 
the hull flow model which, in principle, can treat both the axial 
and crossflow types of separation. In practice, the method was 
implemented in the computer program only in the crossflow-type 
separation case and can handle angles of attack, or sideslip, 
from approximately 30° to 90°. 

Current separation prediction methods can be divided into 
two groups. The methods of the first group require the solution of 
the boundary-layer equations and are either of the differential 
or integral type. The methods of the second group predict sepa- 
ration based on simple formulas or empirical relations and are 
usually very fast and easy to apply. A summary and comparison of 
some of the methods in each group is contained in reference 10. 

One method from the latter group that lends itself particularly 
well to the present situation is the method of Stratford (ref. 

11) which uses only the pressure distribution to predict boundary" 
layer separation. Although originally developed for two-dimen- 
sional flow over a flat plate, it has been shown by Spangler and 
Mendenhall in reference 12 to be effective for bodies in cross- 
flow with only a slight modification. Although the work of ref- 
erences 11 and 12 considers both laminar and turbulent separation, 
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only the latter case is considered for the airship application. 

The Stratford criterion is based upon the concept of divid- 
ing the boundary layer into inner and outer layers. By matching 
the two solutions at the juncture of the layers an equation is 
developed of the following form. 

, dc_ •* x / z .... 


f d Cp 
C P “ "ds 


(10~*Re g ) “V 1 0 = F(s) (14) 

/ 


where s is a streamwise distance parameter. The value of F(s) 
is found by experiment and, as noted by Stratford in reference 
11, is independent of C p . 

The analysis that leads to equation (14) is based on the 
development of a boundary layer on a flat plate with a constant 
pressure imposed up to a distance s^, followed by an adverse 
pressure gradient beyond s m . For purposes of determining the 
parameter F(s), Stratford uses a virtual origin at a position 
s' such that an all-turbulent boundary layer starting at s' 
has the same thickness and velocity profile at s m as the actual 
laminar-transitional-turbulent boundary layer starting at the 
plate leading edge (s=0) . For the fully turbulent flow, refer- 
ence 10 defines the location of the virtual origin (after drop- 
ping a very small term) as: 

S ** 

r m 


r rn (V \ 

- . JT 

TTW 


For the present case of crossflow over a circular body, a 
favorable pressure gradient occurs from the stagnation point to the 
maximum velocity point on the flank. Application of equation 
(15) provides a virtual origin for the circular body case, with 
U m the velocity at the flank of the body. For the circular body, 
a new pressure coefficient has to be defined based on the maximum 
velocity U . 
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( 16 ) 


Using equations (14), (15), and (16), the separation point 
can be predicted as soon as an appropriate value for P(s) at 
separation is selected. Based upon experimental evidence Strat- 
ford found separation to occur at the point where F(s) ■ 0.4. 

The results of reference 12 indicate good agreement with other 
separation prediction methods and data using this value. 

For flow over inclined bodies of revolution at angles of at- 
tack less than 90°, the boundary layer has both an axial and a 
crossflow component. The longer run lengths tend to cause sepa- 
ration closer to the maximum velocity point on the flank of the 
body. As an approximate correction for angle of attack, the fac- 
tor F (s) was replaced with F(s) sin a in the work of reference 
12, with considerable success. That approach is used in the pre- 
sent analysis for angles of attack down to about 30°. Thus, the 
right hand side of equation (14) becomes F(s) sin a. 

In modeling the hull of an airship, the surface is divided 
axially into a number of rings, each of which is divided circum- 
ferentially into a number of panels. The foregoing analysis for 
separation is applied to each ring along the length of the hull to 
locate (usually) two separation points. 

Incorporation of this separation criterion into the potential 
flow model for the hull requires the computer to be able to make 
several nontrivial decisions. One such decision involves an air- 
ship hull in hover. Since there is no free stream, there is no 
well defined stagnation point in the normal sense. However, the 
pressure distribution calculated from purely potential flow con- 
siderations indicates the likelihood of flow separation in the 
vicinity of the rotors. But to apply the separation criterion 
the location of the stagnation point must be specified (i.e., the 


17 


stagnation point is the lower limit in the integral of equation 

(15)). Hence, when applying the separation criterion for a hull 

in the presence of rotors the point s = 0 (stagnation point) is 

defined as the location of maximum C . For the case of one ro- 

P 

tor on a hull in hover, this point will be located on the side 
the hull opposite to the rotor. 

Since the inviscid pressure distributions which govern the 
onset of separation can be very steep for a powered hull in hover, 
a smoothing function is applied before taking the derivatives in 
equation (14) . This is accomplished by fitting a cubic spline 
function to the pressure distribution and then using the coeffi- 
cients of the cubic spline to calculate the derivatives. 

Pressures in Separated Region 

Because of asymmetries in the pressure distribution around 
the hull for all cases except the hull alone out of ground effect, 
the separation criterion is independently applied in both direc- 
tions along each ring starting at the stagnation point. Thus 
two separation locations are determined independently of one an- 
other. The pressures at these two locations are then averaged .. 
give the pressure in the separated region. The pressure on eacn 
panel in the separated region is then given this value. Thus the 
final pressure distribution is the inviscid pressure up to the 
separation points and the average inviscid pressure of the two 
separation points at all points downstream of separation in the wake. 

CALCULATION PROCEDURE 


For ease of computation all necessary computer programs for 
calculating loads on a complete HLA configuration have been com- 
bined into a single operational unit named HELIAIM (Heavy-Lifter 
Aerodynamic Interference Model). This includes a geometry pro- 
gram which calculates the locations of the corner points of the 
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hull source panels, the hull potential flow program, the rotor 
wake program, and the hull flow separation program. In this sec- 
tion the calculation sequence is outlined for a complete airship 
configuration with separation effects included. 

The calculations proceed as shown on the following page. 

The first step is to determine the corner locations of the body 
source panels on the hull surface. A geometry program has been 
developed to generate these coordinates for a body of revolution 
consisting of a cylindrical center section of arbitrary length 
and ellipsoidal end sections. Unequal panel spacing both axially 
and circumferentially may be used if desired. These comer points 
are then used to establish the system of plane quadrilateral 
panels that represent the surface. The no-flow boundary condi- 
tion is sa.isfied at one point on each panel - the control point. 
The panel centroid is taken as the control point. Thus, the next 
step is to determine the velocity induced at each of the control 
points by the rotor system. 

For each rotor, its location, attitude, size, thrust coeffi- 
cient, and angular rate are input, as indicated previously. The 
wake centerline of the rotor is located in the plane perpendicu- 
lar to the rotor disk which contains the onset flow at the rotor 
center. This onset flow could be the free stream only or it 
could be the vector sum of the free stream and a hull-induced 
perturbation velocity. In the calculations made to date, only 
the free stream has been used. The skew angle is calculated to 
locate the wake centerline in this plane. Induced velocities at 
every hull control point are calculated from the rotor singular- 
ity model for each rotor present. 

The rotor- induced velocities are then added to the free- 
stream velocity at each control point and the required source 
density distribution is determined. Once this is done the re- 
sultant flow velocities are calculated by adding the contributions 
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START 


1. Geometry program (GEOM) sets up 
corner locations of body source 
panels on hull surface 


2. Hull program (FATBOD) sets up body 
source panels and calculates location 
of control points 


3. Rotor wake program (HELIWAK) eval- 
uates induced velocity at each control 
point due to the rotor wakes 


4. Hull program takes induced velocity 
plus free-stream velocity at each con- 
trol point and solves for source den- 
sity distribution on body surface 


5. Calculate surface velocities from 
source distribution 


6. Determine inviscid pressure dis- 
tribution from Bernoulli's equation 


7. Evaluate effects of flow separation 
on hull surface with separation pro- 
gram (SEPRaT) 


8. Modify pressure distribution to in- 
clude effects of separation 


9. Calculate loads on hull from modified 
pressure distribution 


Sequence of Program Computations 


2 







of the source elements, the rotors, and the free-stream velocity. 
These velocities are next used to compute the inviscid pressure 
distribution over the hull surface from Bernoulli's equation. 

The inviscid pressure distribution is next evaluated for 
possible flow separation according to the modified Stratford cri- 
terion described in the previous section. Separation effects are 
evaluated around the hull at each axial location, one rinq at a 
time. If separation is found to occ»r, the pressure distribution 
downstream of separation is modified and the program moves on to 
the next ring. This complete modified pressure distribution is 
then used to calculate the total loads on the hull due to effects 
of free stream, rotor wakes, and flow separation. 

At this point one could improve the calculations by iterat- 
ing on the wake location that was determined only from the free- 
stream velocity. The hull program gives the value of the veloc- 
ity at points away from the hull surface induced by the presence 
of the hull. Thus it would be possible to evaluate the velocity 
at the location of the rotors induced by the hull and other rotors, 
re-calculate the wake trajectories and the wake induced velocities 
on the hull surface, and repeat the set of computations until the 
loads converged to a final value. Thi«* procedure has not been 
incorporated in the program at this time. 

The presence of a ground plane can be incorporated by making 
use of symmetry provisions in the Douglas program. The only dif- 
ficulty lies in handling the rotor wake-induced velocities pro- 
perly. In the ground plane case, an image rotor exists for each 
real rotor. The program . ccounts properly for the image hull 
source contribution to the velocities at a control point, so an 
additional feature is used to add on image rotor induced velocity 
as well. Otherwise, the calculation proceeds as before. 
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RESULTS 


In this section, a number of results are shown to indicate 
the capabilities of the method and the types of loading that can 
be induced by rotors and rotor wakes adjacent to an airship hull. 
For purposes of illustration, the analytical methods are applied 
to the "Heavy Lift Airship" (HLA) configuration that evolved from 
the NASA sponsored study reported in reference 1. An additional 
reason for using this configuration is that some small scale mod- 
el data on loads are available from reference 3. 

A drawing of the HLA configuration is shown in figure 4. 

The dimensions indicated are full scale. The hull is circular 
in cross section and has a cylindrical center section and ellip- 
soidal ends. The four rotors are mounted relatively low on the 
cylindrical portion of the hull. The positive x-direction is 
towards the bow, although the configuration is symmetric fore 
and aft. 

The predicted results obtained on this configuration are 
compared with the data of reference 3 wherever possible. These 
are the only data of this kind available, to our knowledge. As 
a general comment on these comparisons, it is noted that consid- 
erable difficulty was experienced with the instrumentation in the 
wind tunnel tests. One overall and four rotor six-component bal- 
ances were used in the tests, generally measuring loads of a few 
pounds or less. The hull loads were obtained as the difference 
between the total load and the sum of the four rotor loads. A 
good bit of data was lost or suspect because of balance failure, 
vibration, and amplifier drift at the low signal levels. As a 
result the accuracy of the data is somewhat questionable, al- 
though the general levels and trends ought to be good. 
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Hull Alone 


For purposes of modeling the hull, its surface was divided 
into panels, with 16 unequally spaced along the length and 24 
equally spaced around the circumference. The panel arrangement 
is shown in figure 5. This arrangement is a compromise between 
a large number of panels to get detailed pressure distributions, 
particularly near separation, and a small number of panels to 
minimize computation times. From an accuracy standpoint, one 
could use fewer panels, particularly on the cylindrical section. 

This is shown in figure 6, which compares pressure distributions 
for the 16 * 24 arrangement with those for a 6 * 8 arrangement 
for the hull in crossflow. The two agree well except on the 
ellipsoidal nose. 

For the hull alone out of ground effect, a potential flow 
solution would indicate no net force for all combinations of 
angles of attack and sideslip. Thus all loads are due to the influ- 
ence of flow separation on the hull. Figure 7 shows the predicted 
separation lines on the hull for zero angle of attack and various 
sideslip angles. At 8 = 90°, the separation line is about 130° 
from the windward meridian except near the very ends of the hull 
and is symmetric fore and aft. At smaller sideslip angles, the 
separation line shifts such that the nose remains unseparated and 
more of the aft end is in separated flow. 

Measured and predicted side force and yawing moment variation 
with sideslip angle for the hull out of ground effect are shown 
in figure 8. The data are from reference 3 and are based on a 
roughened hull which was found necessary to ensure turbulent sep- 
aration. The predicted side force is lower than the measured 
results except at 8 = 90 but shows the same trend as the data. 

The predicted yawing moment also shows the same trend as the data 
but is somewhat low in the 30° to 60° range. 
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The influence of ground effect is shown in the next two 
figures. Figure 9 shows circumferential pressure distributions 
near the center of the hull for the 90° sideslip case with and 
without ground effect. Over the top of the hull, the pressures 
are very similar, and separation occurs about a degree further 
around the hull. Over the lower surface, the ground effect causes 
a higher velocity and a more negative peak pressure coefficient. 
The increased adverse pressure gradient causes separation to occur 
2° to 3° earlier. The average of the two (top and bottom) pres- 
sures at separation in ground effect is a somewhat lower pressure 
coefficient than for the out-of -ground effect case. 

Figure 10 shows the side-force coefficient variation with 
sideslip at H/D=l. In this case the trend and the predicted 
valuer are in reasonable agreement with the data. Because of the 
lower separation pressure the predicted side force increases a 
small amount due to ground effect (compare figures 8 and 10) . 

The data indicate that the influence of the ground on side force 
is small but no trend is evident. 

Rotor Alone 

To validate the accuracy of the rotor wake program some cal- 
culations were made for a single rotor operating in a uniform flow 
and compared with experimental results. This section discusses 
these comparisons. Although a great deal of research, both theo- 
retical and experimental, has gone into the accurate representa- 
tion of the rotor wake, the purpose of most of this research has 
been the accurate prediction of induced inflow at the rotor disk. 
The result of this is that existing induced velocity data is in 
the vicinity of the rotor disk and the wake and not in the rela- 
tively distant field where the HLA hull would be located. Thus, 
the comparisons in this section are not at locations representa- 
tive of the HLA hull. The assumption is made, therefore, that 
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accurate representation of the rotor wake near field implies ac- 
curate representation of the far field. 

The rotor wake analysis assumes that the wake skew angle is 
known a priori. As described in Appendix A the skew angle is 
evaluated from momentum theory for the particular combination of 
rotor geometry, rotational speed, loading, and forward velocity. 

The value derived from momentum theory is compared to experimental 
results taken from reference 13 in figure 11. The agreement is 
quite good, with skew angle being predicted within a few degrees 
by momentum theory. As a point of reference, in the HLA config- 
uration results to be discussed subsequently, for a typical disk 
loading of 478N/m 2 (9 lbs/ft 2 ) in a free-stream velocity of 
15.4m/sec (30 knots) the skew angle is approximately 56°. 

For a given rotor, the total induced velocity at any point 
is the sum of the contributions from the four individual vortex 
systems. Figure 12(a) shows these relative contributions at a 
point outboard and just above the rotor disk (the nearest hull 
control point of the HLA configuration) . Note that the four com- 
ponent velocities in figures 12 are not additive to get the total 
because chey are magnitudes of vectors which have different direc- 
tions. As can be seen, the major contribution is from the vortex 
rings. The difference between the total induced velocity and the 
induced velocity due to the vortex rings is generally within five 
percent of the free-stream velocity for all velocities of interest. 

This difference is even less at points further away from the 
rotor disk. The total induced velocity and that induced by the 
vortex rings at a point near the nose of the HLA configuration 
are compared in figure 12(b). Here the difference is essentially 
negligible, indicating that for points in the far field of the 
rotor wake the wake-induced velocity can be adequately represented 
by accounting tor only the ring vorticity (i.e., using an actuator 
disk model for the rotor wake) . 
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Comparisons of induced velocity are made in figure 13 with 
experimental results from reference 14 for a rotor alone in a 
free stream flow. These curves show the variation in the compon- 
ent of induced velocity normal to the rotor disk, normalized by 
the induced velocity derived from momentum theory, across the 
plane of the rotor disk at two streamwise positions and two heights 
above the disk. Agreement is good at most locations across the 
disk, although the predicted results tend to show more right-to- 
left asymmetry due to advancing/retreating blade aerodynamics 
than do the data. The free-stream velocity for the particular 
set of conditions analyzed in figure 13 is 6.02 times the jet 
velocity. 

Since these experimental results show only the normal com- 
ponent of induced velocity, some calculations were made of induced 
velocities at a number of points in proximity to the rotor wake 
to determine the general nature of the induced flow. Figures 
14(a) and 14(b) show the results at points one-half rotor radius 
away from the rotor wake and at 90° intervals about the wake 
starting at the x-axis. The velocity vectors shown are scaled 
relative to the free-stream velocity. The thrust and flow condi- 
tions used are typical of the rotor used in the HLA configuration. 

In figure 14(a), the side view shows the rotor-induced ve- 
locities directed roughly normal to the wake axis, except in the 
vicinity of the rotor disk. Upstream and downstream of the wake, 
the induced velocities are directed upstream, whereas outboard 
of the flanks of the wake cylinders, the induced velocities are 
directed downstream. This pattern generally resembles that in- 
duced by a doublet which is used in two dimensional potential 
flow to represent the induced velocity of a solid cylinder in a 
uniform flow, as shown in the sketch. 
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Thus, the singularity representation of the wake tends to make 
the wake resemble a solid cylinder around which the free stream 
passes. 

The flow near the rotor disk is generally into the disk, as 
would be expected. Essentially, the rotor disk acts as a sink to 
attract the flow from all directions. 

Figures 14(a) and 14(b) together indicate that the induced 
velocity is higher on the positive y R side of the disk, which 
is the advancing blade side. 

These induced velocity characteristics are useful in analyz- 
ing the aerodynamic interference effects induced on the hull by 
adjacent rotors which are discussed in the following section. 

Hull-Rotor Combinations 

Results for five cases are presented to indicate the nature 
of the predicted interference effects between the rotors and hull 
for the HLA configuration. The five cases consist of a hover 
case out of ground effect, an axial free-stream flow in and out 
of ground effect, and a crossflow case in and out of ground effect. 

The first case is the HLA configuration in zero free stream 
flow with a single rotor on the forward left side developing 
thrust with a disk loading of 395 N/m" (8.26 psf ) . This case 
corresponds to a test condition run in the wind tunnel program 
and also illustrates well the interference effects of a single 
rotor. Since there is no free-stream flow, all hull loads are 
induced by the rotor. 
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The velocities induced by the rotor at the panel control 
points on the hull at an axial station 3m (9 feet) aft of the 
rotor-center axial station are shown in figure 15. The vector 
scale length is shown relative to the velocity \T in the rotor 
plane. As would be expected, the rotor disk and wake act as a 
sink and induce flow towards the rotor. Over most of the hull, 
the induced velocity is less than 0.1 Vj , but adjacent to the 
rotor, the highest induced velocity is about 0.3 . 

Predicted circumferential pressure distributions at three 
axial stations are shown in figure 16. The pressure coefficient 
is nondimensionalized by a dynamic pressure corresponding to a 
unit velocity. The station that corresponds to the induced veloc- 
ities is shown by the circles. Over the top of the hull where 
the velocities are small (45° < 0 < 270°), the pressure coeffi- 
cients are small. As one approaches the high velocity region op- 
posite the rotor (6 - 330°), the pressure coefficients reach very 
large negative values (9 = 310° and 350°). Since the induced flow 
along the hull is towards the 9 = 330° point on both the lower 
and left portions of the hull, a "stagnation point" should occur 
where these flows meet, which is shown by the rise in pressure 
coefficient to about -40 at 0 * 330°. Directly opposite the 
rotor center, the pressure coefficient should rise to zero, but 
at the 15 meter station, there is an axial component of the flow 
which prevents this value being reached. 

At x ~ 32 meters, which is near the forward edge of the 
rotor, the pressure coefficients show an acceleration of the flow 
near the rotor (0 = 320°), but not the "stagnation point" beha- 
vior of the 15 meter station. The flow is more affected by the 
axial component in this case. At x = -5 meters, which is well 
aft of the rotor, the flow is again principally forward axially 
along the hull, and the pressure coefficients are small. 

No flow separation calculations were made for this case. It 
is possible with the adverse pressure gradients that exist from 
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310° to 320° and 350° to 340° at the 15 meter station that separa 
tion would occur over perhaps a 20° band; this would be character 
istic of axial stations near the rotor center (x = 17m) . The 


pressure recovery to C 

= - 40 at 6 = 330° 

would be reduced 

somewhat, if separation 

does occur. 


The integrated hull 

loads for this case 

and the correspond 

ing experimental values 

from reference 3 are 

shown below. 

Quantity 

Predicted 

Measured 

C Y 

-.053 

+ .472 

C z 

+ .029 

-.002 

C m 

-.011 

not reported 

C 

n 

-.022 

not reported 


The predicted values are all consistent with a low pressure 
region on the lower, left, forward part of the hull (adjacent to 
the rotor); that is, a downward normal-force, a left side-force, 
a nose-down pitching-moment and a nose- left yawing-moment. The 
data show signs for the Cy and C opposite those of the pre- 
dicted values. No moment data are available for a consistency 
check. The reason for the differences is not known. The hull- 
alone data are obtained by differencing the total load from the 
main balance and the sum of the four rotor balance loads, and in 
most cases the difference is small compared to the two numbers 
being differenced. Also the hull-alone data include the effect 
of the rotor support arms, and the size of this effect is diffi- 
cult to assess. These may be contributing factors to the 
differences . 

The second and third cases are a hull with four rotors in an 
axial flow (a = 6 = 0) of 15.4m/s (30 knots). In the second 
case, the vehicle is out of ground effect and in the third case, 
the hull centerline is one hull diameter above the ground. Out 
of ground effect, the forward rotors (Nos 1 and 4) have a disk 
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loading of 609 N/m 2 (12.73 psf) and the aft rotors 521 N/m 2 (10.88 
psf ) , which are the values measured in the wind tunnel for this 
case. The resulting wake skew angles are 51° and 54°, respectively. 

In ground effect, the disk loading of all four rotors is 
479 N/m 2 (10 psf) , again in consonance with the tunnel tests, and 
the skew angle is 56°. 

The wake configurations are shown in figure 17 for wake skew 
angles of 56°. In the second case (without ground effect), there 
is no image system and the rotor wakes are swept slightly forward 
of those shown in the figure. In the third case, the image sys- 
tem exists, as shown. The sketch shows a limitation in the use 
of the Heyson rotor model, which because of the analytical inte- 
grations performed, is constrained to a straight centerline. The 
image wakes from the forward rotors come close to the aft end of the 
hull and tend to cause unrealistic induced loads. A more realis- 
tic wake model would be curved to approach and become parallel to 
the ground. To go to a curved wake would require dropping the 
analytical integration and modeling the wake with a series of vor- 
tex rings and axial finite-length vortex filaments, with numerical 
summation of the induced velocities, as was done in reference 15. 

Rotor-induced velocities in and out of ground effect are 
shown in figure 18 for an axial station near the center of the 
hull between the fore and aft set of rotors. There is some in- 
fluence of the image system on the velocities, but the two sets 
are generally similar. The rotors induce a downwash over the top 
of the hull, an upwash over the bottom of the hull, and outward 
flow towards the rotor wakes near the rotor disks. The flow is 
not symmetrical right-to-left because the rotors are all rotating 
in the same direction and their induced velocities in a crossflow 
are different on the advancing and retreating sides. 

The corresponding circumferential pressure distributions are 
shown in figure 19. Over the top of the hull, (0° < 0 < 180°), 
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the pressure is relatively uniform and not greatly different from 
the hulL-alone value. There are small pressure coefficient peaks 
(least negative) at 0 « 210° and 330° which correspond to the 
crossflow velocity reversal along the hull opposite the rotors. 
Over the bottom portion of the hull there is a suction peak in- 
duced by the high rotor-induced velocities. 

The integrated hull loads for these cases and the measured 
values from reference 3 are shown below. No predicted flow sepa- 
ration effects are included in these calculations. 

Predicted Values Measured Values 

Quantity Out of G.E. In G.E, Out of G.E. In G.E. 


-.006 

+ .088 

not reported 

not 

reported 

+ .156 

+ .272 

-.443 

- 

-.035 

-.024 

+ .068 

not reported 

not 

reported 

-.135 

- .097 

not reported 

not 

reported 


Out of ground effect, the predicted centers of pressure for 
both the vertical and lateral forces are forward of the hull cen- 
ter, as was the case for the single rotor-hull in hover. In the 
presence of the ground, the side force changes sign and its cen- 
ter of pressure moves aft of the hull center, since the sign of 
the yawing moment is not changed. The vertical force maintains 
the same sign, increases in magnitude and its center of pressure 
moves aft of the hull center, since the sign of the pitching mo- 
ment changes. The aftward shift of the centers of pressure is 
probably due to the influence of the image wake on the flow over 
the aft end of the hull. 

Only the normal-force results are available from the tests. 
Although the change in force due to the ground presence has the 
same sign as the prediction, the data indicate an opposite sign 
to the predicted values. Again the reasons for the differences 
are not known . 
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The fourth and fifth cases are a hull with four rotors in 
a crossflow (a ■ 0, 0 « -90°) of 15.4 m/s (30 knots). The fourth 
case is for the vehicle out of ground effect and the fifth is for 
the hull centerline one diameter above the ground. All four ro- 
tors have a disk loading of 478 N/m 2 (10 psf) out of ground effect 
and in ground effect; the wake skew angles are 56°. 

The wake configuration is shown in figure 20. In this case, 
the image wakes do not intersect the hull, but they come fairly 
close and do induce some large velocities on tha downstream sur- 
face of the hull. As indicated previously, a curved wake would 
provide a more accurate simulation. 

Rotor induced velocities at an axial station near the aft 
rotor center are shown in figure 21. Without ground effect, the 
rotors induce a downward flow over the hull which is largest near 
the two rotor disks. In the presence of the ground plane, the 
image wake induces some large velocities directed away from the 
wake generally towards the free stream flow. Thus, large differ- 
ences in induced velocities occur due to the ground plane effect. 

The resulting pressure distributions are shown in figure 22. 
There are three curves shown: a hull alone out of ground effect 

(without separation) and the hull with rotors in and out of 
ground effect (with separation) . Over the top of the hull 
(0 < 0 < 120°), the three curves are similar. Separation is pre- 
dicted to occur at about 130° for both in and out of ground effect. 
The pressure distribution over the lower side of the hull up to 
separation is considerably different than the hull alone. The 
influence of the rotors is to induce a flow counter to the free 
stream flow over the lower part of the hull, which yields a high- 
er pressure than for the hull alone. Separation occurs on the 
lower part of the hull about 130° around from the windward merid- 
ian for both the in and out of ground effect cases. 
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It is of interest to note in connection with figures 21 and 
22 that while the rotor-induced velocities show the reversal in 
direction along the hull that was noted in connection with the 
hover case (figures 15 and 16) , when the free stream is added, the 
reversal in total velocity along the surface does not occur. 

That is, the pressure distribution with a free stream (figure 22) 
does not show the "stagnation point" behavior that the hover case 
shows (figure 16) . As the disk loading increases and/or the free 
stream decreases, the pressure distribution would proceed from 
the more regular type of figure 22 to the more complicated one of 
figure 16, with or without ground effect. 

The predicted and measured loads on the hull for these condi- 


tions are 

shown in the table below. 

A large side 

(or drag) force 


Predicted 

Values 

Measured Values 

Quantity 

Out of G.E, 

In G.E. 

Out of G.E. 

In G.E. 

C Y 

1.59 

.74 

2.22 

2.77 

C 

z 

-.18 

-.26 

1.09 

1.54 

C m 

-.067 

-.083 

not reported 

not reported 

C 

n 

.05 

+ .57 

not reported 

not reported 


is predicted due to separation out of ground effect which is 
decreased in the presence of the ground. This effect is due 
to the proximity of the image wake and would probably be much 
less significant with a curved wake model. The measured side 
force out of ground effect is somewhat higher than the predicted 
value. Since the disk loading changed somewhat for the tests in 
ground effect, the force changes due to presence of the ground 
are not comparable to the predicted changes. The predicted normal 
forces are small and upward in direction, because of the absence 
of the suction peak on the lower hull surface. Again there is a 
difference in sign between predicted and measured normal force, 

which cannot be explained. 

* 

These data correspond to a disk loading of 539 N/m (11.25 psf) . 
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Systematic Rotor Effects 

One purpose of developing this analytical method to predict 
the effects of aerodynamic interference on airship configurations 
is to use the method to investigate various design changes, such 
as hull shape or rotor placement. By studying the results of 
such changes not only can an optimum design with least unfavorable 
interference be determined but much can also be learned about the 
basic interaction mechanisms and how they affect the total vehicle. 

A few results on the effects of rotor placement are presented 
in this section. The baseline HLA hull (figure 4) was used 
throughout and only the vertical and lateral locations of the 
rotors were changed. The axi^l locations remain the same as in 
figure 4. These results are for a vehicle out of ground effect 
in a crossflow (a » 0° , $ ■ -90°) of 9.14 m/sec (30 ft/sec), with 
a disk loading of 478 N/m 2 (10 psf), and with separation on the hull. 

The effect on hull loads of vertical placement of the rotors 
is shown in figure 23. The horizontal distance between the tip 
of the rotor and the surface of the hull has been kept constant 
at the same value as the baseline configuration in figure 4. 

Side force and yawing and pitching moments are only slightly af- 
fected by rotor vertical placement. The normal force, however, 

is strongly affected, changing from an upward force (negative C ) 

z 

for high rotor positions to a downward force when the rotor is 
even with the bottom of the hull (Z Q /R = 1) . Examination of the 
induced velocity fields revealed that when Zq/R = 1 the rotor- 
induced velocity over the hull was directed downwards everywhere 
and was very large in magnitude. Conversely, when Zq/R = 0 
there was a region of upwash on the lower part of the hull and 
the magnitude of the downwash was much less. Hence, the vertical 
rotor vertical placement has a substantial effect on the rotor- 
induced velocity on the hull and, subsequently, on the hull ver- 
tical loads. For this example, a change in C from -0.6 to 0.4 

z 

with rotor height change corresponds to a change in hull vertical 
force approximately equal to the thrust of one rotor. 
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Effects due to changes in lateral spacing are shown in figure 
24. In this study the vertical location of the rotors is 
Zq/R ■ 0.5. Normal force, yawing moment, and pitching moment all 
decrease as the rotors are moved farther away, while the side 
force overshoots and then decreases towards the hull-alone value. 
Thus, the hull loads are also strongly affected by the lateral 
location of the rotors. 

The results in this section, although limited, indicate the 
importance of systematic studies of rotor locations, under many 
different flow conditions, to properly locate the rotors relative 
to the hull to either minimize adverse rotor interference or ob- 
tain a favorable influence. 

CONCLUDING REMARKS 


An analytical method has been developed to predict the per- 
turbation pressure distribution induced on an airship hull by 
adjacent rotors. The method accounts for flow separation on the 
hull at angles of attack and sideslip above about 30° , where there 
can be considerable interaction between the rotor induced pertur- 
bation effects and the separation-induced hull loads. A number 
of cases were calculated and are reported herein to illustrate 
the nature of the hull-rotor interaction. This section contains 
some comments on the work and results. 

One important limitation of the method is the constraint im- 
posed by the Heyson analytic wake model of a straight, semi-infi- 
nite wake. Out of ground effect, the limitation xs probably of 
little importance. The significant wake-induced effects are 
caused by that portion of the wake in close proximity to the hull. 
It is relatively simple to modify the wake skew angle to best lo- 
cate the wake relative to the hull considering both the free stream 
and the hull flow perturbation. This approach should suffice for 
all configurations except those in which the rotois are very close 
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to and high on the hull and have wakes with many diameters of 
length close to the hull. However, in ground effect, it is felt 
that a straight wake is inadequate, and a curved wake model is 
necessary. In this case, one must forego the analytic model and 
use a numberical model such as that of reference 15. This 
numerical model uses vortex rings concentric with a centerline 
which may be curved arbitrarily. Since the results of figure 12 
show the vortex rings of the Heyson model to be the major con- 
tributor t ' rotor-induced velocities on the hull, the numerical 
model shou’ 1 be satisfactory. The vortex ring model has the 
additional feature of modeling entrainment through change in the 
diameters of the rings along the wake, should that be desirable. 
For high crossflow, a surface vortex quadralateral distribution 
can be added to model blockage, which will show much the same 
characteristics as in figure 14. This model has been implemented 
in other work at Nielsen Engineering and can be readily incorpo- 
rated in the airship-rotor interaction model. 

Little was said of the influence of the hull on the rotor 
performance. In the cases examined, the experimental thrust coef- 
ficients from reference 3 were used. These tests showed only a 
small influence of hull presence on rotor thrust (generally less 
than 10 percent) , even for the case of a rotor operating down- 
stream of the hull in crossflow. There was a greater influence 
of ground presence on rotor thrust (as high as 20 percent) . Some 
effort needs to be made to estimate these effects and develop a 
systematic approach to calculate the change in rotor thrust as 
the rotor is brought into proximity with the hull, in and out of 
ground effect. 

Due to limitations in scope of the effort, no calculations 
were made to illustrate the capability of the method to model a 
hull with tail surfaces to examine the influence of tails on hull- 
rotor interference. The Douglas program has been used by NEAR to 
make calculations for deltoid airship configurations with tails. 
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although no rotor-induced effects were included. The method 
would be limited to angles of attack and sideslip for which no 
separation occurs on the tails, although one would probably not 
operate an airship with tails at very large angles, at least in 
steady state. 

The approach that was used to calculate the pressures on 
the hull in l£»rge regions of separated flow is rather straight- 
forward. The approach works well for the hull alone over a large 
angle range and is felt to be a good initial approach to model 
this phenomena, which can be important in many cases. The method 
becomes less satisfactory when the two pressures at separation on 
a given ring begin to differ substantially due to rotor-induced 
effects, since these two pressures are averaged. In order to get 
a better understanding of this problem and develop a more accurate 
model, it is essential to perform experiments in which pressure 
distributions are measured on a hull with rotor interaction over 
a range of disk loadings, free stream velocities, and angles of 
sideslip. 

The discrepancies in sign between predicted and measured 
normal force need to be resolved. One of the questions that arose 
in the work was whether or not the hull source singularity model 
was providing a zero net source strength in cases of nonuniform on 
set flow. For the hull-alone out of ground effect, one sees a 
symmetry fore-and-aft or circumferentially, depending on the direc 
tion of the onset flow, which provides zero net source strength. 
Since the Douglas program does not require zero net strength as a 
condition on calculation of source strengths, the computer program 
was modified to sum the source strengths after they are calculated 
to investigate this question. Results show that with a nonuniform 
onset flow, the surface source strengths are redistributed non- 
uniformly to cancel locally high onset flow velocities in such 
a way that the zero net source strength condition is preserved 
over the entire body, provided it is closed. Since no other 
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fundamental issues arise in application of the methods to this 
problem, the differences in predicted and measured normal force 
are attributable either to the separation and wake pressure 
aspects of the calculation or to problems with the small scale 
data. Only carefully executed experiments to obtain surface 
pressure distributions and surface flow visualization will resolve 
this problem. 

Finally, the Douglas program is not well suited in its ori- 
ginal form to the hull-rotor problem because its emphasis is on 
single body problems (as opposed to ground-plane image cases) 
with predominately lifting-surface configurations. Program 
modifications have been made to handle every rotor-airship case 
of interest, but many cases are inefficient to run and it is 
very difficult to prepare input. Additional work needs to be 
done to simplify the program and make it much easier to prepare 
input and interpret the output. 

In summary, the recommendations for additional investiga- 
tions include the following: 

• experiments at high Reynolds numbers on typical hull-rotor 
configurations in which simultaneous measurements of rotor 
and hull loads, surface pressures, and surface flow 
visualization are taken 

• implementation of a curved wake rotor model to improve 
the prediction method for cases in ground effect and with 
closely coupled rotors. A part of this task is the defi- 
nition of the wake centerline location for flow conditions 
of interest 

• development of methods for estimating the influence of 
the hull on thrust of rotors adjacent to the hull 

• improvement in the separation location and wake pressure 
models to include nonuniformities in the wake flow region, 
based on critical experimental results 
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• modification of the computer program to make its operation 
more efficient and to make it user oriented in terms of 
input preparation and output interpretation. 

Nielsen Engineering and Research, Inc. 

Mountain View, California 
April 1978 
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APPENDIX A 


DERIVATION OF ROTOR WAKE MODEL 

As discussed in the main report, the mathematical model used 
to describe the rotor wake is based on the classical vortex theory 
of a non-contracting wake. The model, shown in figure 2 of the 
main report, divides the wake into four vortex systems. The ori- 
gin of each of these vortex systems has been described in the 
main report and will not be repeated here. In this section the 
induced velocity at a point due to each of these vortex systems 
is derived. The coordinate system used in this development is 
given in figure A.l. The onset velocity, by definition, lies in 
the X R , Z R plane. 


FORM OF BLADE CIRCULATION 


As explained in the main report, the blade circulation varies 
in the azimuthal plane according to the relation 

T = I' + I" sin i/i (Al) 

0 1 

Now, the incremental thrust along the blade is 

dT = U T r dr (A2) 

where U T is the tangential component of the resultant velocity 
at the blade element, i.e. 


U T = i*r + u ilR sin ij’ 


(A3) 


Therefore, the rotor thrust is given by 


T = j- 
2tt 


2n 


+ iu2R sin ij 1 ) ( I' ^ + I’ ( sin ip ) drdC‘ 


(A4) 


o o 
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(A5) 


T - bp.Sgi (r 0 ♦ pr,» 

Expressed in terms of the non-dimensional thrust coefficient equa- 
tion (A5) becomes 


C 


T 


T 

pirR 2 (flR) 2 


b 

2 nuR 2 


(A6) 


Defining the thrust moment as 


M,p = I P a> (r<r + pftR sin ifi) (T q + Tj sin iji)r dr, 
and carrying out the integration gives 

3 


(A7) 


Hp - p„^R' 


r 

r r, 

ur 

[T + 

-1 + — 2 . 
[3 2 J 

sin tp + — j— sin 2 ip 


(A8) 


If the assumption is made that flapping terms above the first 
harmonic are negligible, then the thrust moment can be shown (ref- 
erence A . 1 ) to be independent of the azimuth angle Thus, 

from equation (A8) , 



(A9 ) 


or 



(A10) 


Substituting this expression into equation (A6) the following re 
lations for blade circulation in terms of thrust coefficient are 
obtained . 


r 


0 


2C T fiTtR 2 



(All) 
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auc^nR 2 



(A12) 


EVALUATION OF WAKE FKEW ANGLE 


The wake skew angle is defined as 


X 


-l 

* tan 


I±L 

A 


where 


A = 


sin a + V . 

_2 1 

fiR 


P = 


V COS ^ 

00 QJ 


-w 


(A13) 


(A14) 


Hence, the normal component of the induced velocity at the rotor 
center, Vj , is required to determine the skew angle. Using clas 
sical momentum theory to calculate (e.g. reference A.l) the 

following relation is obtained 


-1/2 C AR 
V. * — — - - / 

Vp 2 + a 2 


(A15) 


Substituting the values for y and A , the value of is 

given by 


V 4 + (2V sin a) 


V? + (V 2 ) V 2 + 
j 00 : 



4 ") 

(flR) 

4 


= 0 (A16) 


Equations (A13) , (A14) , and (A16) are thus combined to yield a 
value for the wake skew angle. 
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VORTEX INTEGRAL RELATIONS 


Central or Root Vortex 

This vortex lies along the wake axis and is the superposition 
of the root vortices of all blades. Hence, its circulation strength 
is the sum of the blade bound vortex strengths, i.e., 

r R » (A17) 


To determine the induced velocity of this single vortex filament 
the following form of the Biot-Savart Law is used. 


V r) 



ds * (r-s) 

i ■+ 1 3 

r-s 


(A18) 


From figure A. 2 the following quantities are determined by 
inspection. 

s = (L sin x * ~ L cos x k) (A19) 

ds = (sin x i - cos x k)dL (A20) 

r-s * (x R - L sin x)i + y R 3 + (z R - L cos x)k (A21) 


Substituting equations (A19) , (A20) , and (A21) into equation (A18) 
gives the following three integral equations for the three com- 
ponents of velocity induced by the root vortex. 


U R = 


-br c 00 

0 


4tt 


y„ cos x dL 


o [ (x R - L sin x> 2 + y R + 


(z„ - L cos x) 2 1 3/ ^ 2 


(A22 ) 


V 


R 


br 

o 

4 IT 


sin x(z R “ L cos x) + cos X( X R ~ L sin X) 
o [ (x R - L sin x) 2 + y R + (z R * L cos x) 2 ]* /2 


dL 

(A23) 
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w r 


-br 
< 

4ir 


r 


y R sin x dL 


T7T 


[ (x R - L sin X) 2 + y R + (* R " L cos x) 2 ] ^ 


(A24) 


Outer Wake 

Since the outer wake is composed of both circumferential and 
axial vorticity, the combined strength being the strength of the 
tip vortex, the relation between these two systems must first be 
determined. This is done in the following manner. 

Assume that one vortex element is formed per revolution per 
blade. Then the number of elements formed per unit time is 
bfl/2ir. These elements are carried off with the velocity 

flR V y i + P (A2 5) 

the wake velocity at the rotor. Therefore, the vorticity on 
the edge of the wake is 

-0 *^== ■ £ ■ £ + £r «*» ♦ <A26 > 

ftR Vu 2 + X 2 

In one revolution the circumferential distance traveled by 
a blade tip is 2 ttR and the axial distance traveled is L. The 
relation between the axial and the circumferential vorticity is 
best seen by unrolling one revolution worth of vorticity shed 
by one blade as in the following sketch 



I 

2itR 
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The axial vorticity is <i> a - w ToT sin $, and the circumferential 
vorticity is u> c ■ w ToT cos 3* Hence, the total vorticity is 

i“ ToT i - v»i + - “tot < m7 > 

Vectorially, this is written + u> c ** u> To ,p. Since one revolu- 
tion occurs in 2ir/ft sec., then using the wake velocity at the 
rotor given by equation (A25) the distance L becomes 

L = 2irR V )J 2 + X 2 (A28) 

If L = 0 (ft -► ®) then all the vorticity is circumferential and 
if L •* « (ft -*• ®) then all the vorticity is axial. Thus the 
assumption is made that the ratio of axial vorticity to circum- 
ferential vorticity is in the same ratio as L/2 ttR, i.e., 

!r - SH - (A29> 

C 


or 


u a * u) c ^ p 2 + X 2 (A30) 

Substituting equation (A30) into equation (A27) . 


to 


c 


^ToT 

[1 + U 2 + A 2 ] 7 ' 


(A31) 


Using the expression for the total vorticity given by equation 
(A26) , the circumferential vorticity is written in the following 
form. 


f d rcl 

[■arj 


bft 1 

2 TT f 

ftRVU 2 + A 2 


[1+ p 2 + A 2 ] ^ 


< r „ + r , 

1/2 0 1 


sin i|») 

(A32) 
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The relation for the axial vorticity is found by combining equa- 
tions (A30) and (A32) . 


f dr 1 

bn 

(T + r sin 1|>) 

o . 1 _ ij 2 j. \ 2 

r dr i 

c 


2 TT 

flR[l + y 2 + X 2 ] l/ 

dL 


(A33) 


Having found expressions for the axial and circumferential vorti- 
city, the analysis can proceed to the determination of the induced 
velocities due to these two vortex systems. The induced velocity 
due to the circumferential vorticity will be treated first. 

Circumferential vorticity .- From the Biot-Savart Law, the 
induced velocity can be written as 


dV 


^1 

4 TT 



ds x (r-s) 


I *+■ -> . 3 

r-s 


dL 


(A34) 


From figure A. 3, the following relations can be obtained. 


s = i(R cos ii» + L sin x) + 3 (R sin ij;) + &(-L cos x) 
ds = [i(-R sin <J>) + } (R cos ) ] dip 
r-s * i(x R - R cos \p - L sin x) + 3 (y n “ R sin ty) + £(z. 


(A35) 

(A36) 

+ L cos x) 
(A37) 


Substituting the above expressions into equation (A34), the follow- 
ing relations are obtained for the velocities induced by the cir- 
cumferential vorticity. 

r dr c ^ (z R +Lcosx) cosij; dLdip 

" dL [ ( x R" Rcos ^’“ Lsin X) 2+ (y R -Rsin-4)) 2 + (z R +Lcosx) 2 1 3//2 

(A38) 


U = 
c 4 tt 


> 2 TT f oo 
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-R f 2, (" f dr c) 

^ Jo Jo W 


(z r +Lcosx) sinij; dLdij; 


T7T 


[ (x R -Rcos4»-Lsinx) 2 + (y R -Rsint|>) 2 + (z R +Lcosx ) 2 ] ^ 


(A39 ) 


if n: 


w. 


(x R cos^-y R sin^+R+Lsinxcosij/) 


T7T 

' o ■'o v — ' [ (x R -Rcos^-Lsinx) 2 + (y R -Rsin^) 2 + (z R +Lcosx) 2 ] ' 

(A40) 

Axial vorticity .- For the axial vorticity, referring to fig- 
ure A . 4 , the relevant quantities are found to be 


s = 


ds 


r-s = 


x (L sin x + R cos *10 + 3 (R sin ij>) + k(-L cos x) (A41) 

[l sin x “ it cos x]dL (A42) 

i(x R - L sin \ - R cos + 3 (y R - R sin <Jj) + k(z R + L cos \) 

(A43) 


The form of the Biot-Savart Law for the axial vorticity is 


dV = -jr- 
4tt 


dr. 


dR^ 


ds * (r-s) 


r-s 


d(RiJ>) 


(A44 ) 


Substituting equations (A42), (A43) , and (A33) into equation (A44), 

yields the equations for the induced velocities due to the axial 
vorticity as follows 


m = R Vu 2 + A 

a 


1 + ^2 

4t 


(2 tt 

> 00 

f dr c) 

Jo J 

0 

rar 


{y 0 cosx~Rsin^cosx } dLdiJi 


3 /2 


[ (x D -Lsinx-Rcosi^) 2 + (y R -RsiniJ/) 2 + (z_+Lcosx) 2 ] ^ 


(A45) 
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ITS?] 


4w 


{-x_cosx + Rcos^cosx~z 0 sinx} dLd^ 


'0 '0 


[ (x R -Lsinx-RcosiM 2 + (y R -Rsin\M 2 + (z R +Lcosx) 2 ] 


77 * 7 ? 


(A46) 


■. * ^ n:r 


0 • 0 


c 

dL 


{y R 3inx~Rsin4»sinx} dLdip 


[ (x R -Lsinx-Rcos^) 2 + (y R -Rsini|i) 2 + (z R +Lcosx ) 2 ] 

(A47) 


2 i V 


3/2 


Inner Wake 

The radial vorticity is distributed across the interior of 
the wake, as depicted in figure A. 5. The strength of the radial 
vorticity is the derivative of the strength of the blade vorticity. 


_d_ dT 
dii» I dL 



COS 


(A48) 


Again, starting with the appropriate form of the Biot-Savart Law, 


_ 1 d 

dV r " 47 3* 


fdr)ds * (r-s) 


dL 


dLdt|> 


r-s 


(A49) 


the following re . 1 a\ ions are found from figure A. 5. 

s ■ i(r cos + L sin x) + 3 *r sin + k(-L cos x) (A50) 

ds = (I cos + 3 sin ^)dr (A51) 

r-s = i(x R - r cos - L sin x) + 3 (y R - r sin iM + lc(z R + L cos x) 

(A52 ) 


Substituting these relations into equation (A49) the following 
expressions for the velocity components induced by the radial 
vorticit> . 
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u. 


_1. 

4it 



The 
in terms 
The next 


For 
are four 

1 ) 

2 ) 

3) 

4) 


dr {z R simJ>+Lsinijicosx} dLdrdiJ/ 

dL [ (x R -rcos^-Lsinx ) 2 + (y R -rsin^) 2 + (z R +Lcosx ) 2 ] 3/2 

(A53) 

dr {-z r cos^-Lcosi{j} dLdrd^ 

dL [ (x R -rcos^-Lsinx) 2 + (y R -rsini|)) 2 + (z R +Lcosx) 2 ] 3/2 

(A54) 

dp {y R cos^-x R sinii)+Lsinxslmi»} dLdrd4/ 

dL [ (x R -rcosi|»-Lsinx) 2 + (y R -rsin^) 2 + ( 2 R +Lcosx ) 2 ] 3//2 

(A55) 

above derivations express each of the velocity components 
of integral relations for each of the four vortex systems, 
step is to integrate each of these integral relations. 

INTEGRATION OF THE VORTEX EQUATIONS 

each of the three components of induced velocity there 
terms, due to 

Root vorticity; U R , V R , W R 

Circumferential vorticity; U_, V , W 

c c c 

Axial vorticity; U,, V. U? 

a a a 

Radial vorticity; U r , V r , 


nr* 

; o ; o J o 


m oo 
0 


_d_ 
d ¥ 


m oo 
D 


d 

3 ? 


Each of the total of twelve terms has been expressed above in 
integral form involving single, double, or even triple integrals. 
In this section these terms are reduced to integrals over $> , the 
azimuth angle, or, for the case of the root vorticity which is 
not a function of \p , integrated completely. The remaining in- 
tegrals over are then solved numerically. 
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Root Vorticity 


Integrating first the root vorticity terms, equation (A22) 
can be written 


U 


R 


~ br 0 yR cos x 

4tt 




dL 

[ (x R -Lsinx) 2 +y R + ( z r -Lcosx) 2 ] 3/2 


(A56) 


Expanding the denominator of this integrand and expressing it in 
the form Rj = a + bL + cL 2 reduces equation (A56) to the gen- 
eral form 


r 


dL 


i 



The solution to this integral can be found in reference A. 2 
(integral 2.264, no. 5) and is 


dL = 2 (2cL + b) 


(A57) 


where 


A = 4ac - b 2 


(A58) 


Applying this solution to equation (A56) and evaluating the result 
at the limits yields 


U R - 


_br o y R COSX 


1 + (x„sinx+z 0 cosx)/R r 


4 7T 


[R|-x 2 sin 2 x-z 2 cos 2 x-x R z R sinxcosx] 


(A59 ) 


where = x R + y R + z R and is the distance from the origin of 

coordinates (the center of the rotor) to the field point P(x R , 

y R' Z P. ) • 
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The relation for 


V R is given by equation (A23) and is 


V 


R 



■® -sinx (z R -Lcosx ) -cosx (y R -Lsinx) 
o { (x R -Lsinx) 2 t.Vp+ (z R -Lcosx) 2 1 % * 


(A60 ) 


The denominator is the same as in equation (A56) , Expanding the 
numerator shows there are three terms, the first two independent 
of L and the third term varies with L to the first power. 
Hence, the first two terms can make use of the solution for U R 
and the third term is of the general form 


LdL 



The solution to this integral is given by integral 2.264, no. 6 
of reference A. 2; 


LdL _ -2 (2a -t- bL) 

V*I 


(A61) 


Applying equations (A57) and (A61) to equation (A60) and evalua- 
ting the result at the limits yields the following relation for 


-br 

V R = TrT 


(-z R sin x ~ x p COS x) 


1 + x R sin x + z R cos x)/Rj 
Ta74) 


+ (2 sir. x cos \) 


R p + x R sin x + z R cos x 
__ 


(A62 ) 


The relation for W R , given by equation (A24), is 


-bl' 




y R sinxdL 


R 4lT o ( (x R -Lsinx) 2 +y R + (z R “Lcosx) 2 ] 3//2 


(A63) 


51 


The general form for this integral is the same as the general form 
for U R . Thus the same general solution, given by equation (A57) , 
can be used. Writing the final result directly, 


W R “ 


4tt 


y sin x< 

K 


1 + (x R sinx+z R cosx)/Rp 
[R^-x^sin 2 x-ZRCos 2 x-x R z R sinxcosx] 


(A64) 


Circumferential Vorticity 

The x-component of the induced velocity due to the circum- 
ferential vorticity, given by equation (A38) is repeated here. 


u = =5 

c 4 tt 


f2ir 

r 00 

r dr i 



c 

■ o • 

f o 

dL 

^ / 


(z r +Lcosx) cosip dLd4> 


[ (x R -Rcos4>-Lsinx) 2 + (y R -RsimJ>) 2 + (z R +Lcosx) 2 ] 


./2 


(A65) 


The denominator, grouped into powers of L, becomes 

R 2 + 2(z R cos x “ x R sin x + R cos 4* sin x)L + L 2 (A66) 


where 


R 2 = Rp + R 2 “ 2r(x r cos ^ + y R sin 40 . 

This quantity is the distance between the point on the rotor disc 

at the angle 4) and the field point P(x R , y R » z R ) as shown in 

figure A.l. Recalling from equation (A32) that (dr /dL) is in- 

c 

dependent of L, equation (A65) can be expressed in two terms 
as follows. 


U c * 


-R 

4tt 


r 2 7T 

f dr c^ 


0 

( dLj 

Z_COS Ip 

R J, 


dL 


o [R*+2 (z R cosx-x R sinx+Rcos4>sinx) L+L^ ] 


TJT 


dij; 


R 

p 2 TT 

rdr > 

c 

j cos x cos 4 1 

. 00 

f LdL 

4u _ 

0 

dLj 

1 3 2 

o lR 2 +2 (z_cosx~XpSinx+Rcos4)sinx) L+L 2 ] 




(A67 ) 
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Integrating over L by applying equations (A57) and (A61) the 
expression for U c becomes 


u » zR 

•2^ 

f dr c' 

cosiji (R cosx+z D )diJ) 

c K (Ti&rw 

C 4 TT 

! o 

l dL j 

1 InDO; 

[R +z_cosx~x R sinx + Rcos^sinx ] R 


The remaining two velocity components due to the circumfer- 
ential vorticity (V and W ) are of the same general form as 

c c 

U c . That is the integrals can be broken into two terms, one 
corresponding to equation (A57) and the other to equation (A61) . 
The results of the integrations can thus be written directly. 


V = 

C 4 7T 


^ 2 TT 

r dr a 


C 

'0 

L dL J 


sin 4* (z^+R c cos'x)d4> 


[R +z_cosx-x D sinx+Rcos4'sinx] R 

C K K C 


(A69 ) 


W = 

C 4 7T 


2tt 


dr^'j R+R^sinxcos^-XpCOSv-ypSini^d^ 


dL 


[R +z n cosx-x„sinx+Rcos4>sinx] R 


(A70) 


Axial Vorticity 


Velocity terms for the axial vorticity involve the same 
integrals as for the root and circumferential vorticity, equations 
(A57) and (A61) . Hence, the final expressions for these forms 
are written without derivation. From equation (A45) , 


U = 
a 


rV y 2 + A 2 

4n 


r 2 tt 

r dr i 

[ 

c 

Jo 

dL 

V / 


cosx (y R -RsiniM dip 
[R +sinx (Rcos4)-x_ ) +z D cosx ] R 

C K K C 


(A71 ) 


From equation (A46) , 


v a " 


>v — 


4 TT 


' 2 tt 

f dr > 
[ c 


0 

l ^ j 

[R, 


cosx (Rcos’^-Xp) -z R sinx 


(A72) 


[R c +sin\ ( Rcos4j-x r ) +z R cosx ] R c 
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From equation (A47) , 


W_ 


rV 2 


+ X' 


r2:r 


4 IT 


dL 


sinx (y R -Rsini|>) 

[R c +sinx (Rcosii’-x R )+z R cosx]R c 


(A73) 


Radial Vorticity 

These terms, given by equations (A53) , (A54) , and (A55) in- 
volve integrals over both L and r. Integration over L involves 
the same general forms as the previous integrals over L. The 
integration with respect to r is somewhat more involved. The 
general forms of this integration is given in reference A. 3 and 
is not repeated here. The final forms are, from equations (A53) , 
(A54) , and (A55) , 
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where 


A = x R cos ip + y R sin p (A77) 

B = z R cos x “ x R sin x (A78) 

C = (y R cos ip - x R sin ip) cos x “ z R sin x sin 4* (A79) 

As a final comment, some of the expressions derived above 
involve the number of blades, b. However, when the relations for 
the vortex strength are substituted into these relations the blade 
number drops out and the final expressions are independent of 
blade number, as discussed in the main text. 
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Figure A. 4.- Outer Wake Axial Vorticity Nomenclature 
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(a) Outer vortex wake. (b) Circumferential 

component. 


(c) Axial component. 


Figure 3.- Decomposition of cylindrical wake into circumferential 
(vortex ring) and axial components (from reference 7). 
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, Degrees 

Figure 9.- Predicted circumferential pressure distribution 
near hull center (x ■ -5m) for HLA hull alone at 


-90 


in- and out-of -ground effect. 
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P, Degrees 


Figure 10.- Comparison of predicted and measured side force 
on HLA hull alone in ground effect. H/D 3 1. 




oo'' j at hover 


Figure 11.- Comparison of predicted and measured 
skew angles for a rotor alone. 
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(a) Point in near field 
(0. 25 R, -1. 33 R, 0. 3 R). 


Figure 12.- Magnitude of induced velocity contributions 
at a point outside the wake for vortex systems 
comprising a rotor alone in a uniform flow. 





Figure 13*- Comparison of measured and predicted 
vertical- induced velocity for a rotor alone 
in a uniform flow at V^/V^ - 0. 167. 
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(a) Side view. 


Figure 14.- Induced velocities adjacent to a 
rotor in a uniform flow. 




Rotor 


Figure 15.- Rotor- induced velocities on HIA hull in hover 
at x * 15 m (48 ft). T /S 0 ■ 395 N/m 2 (8.26 psf). 
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Figure 16.- Predicted circumferential pressure distributions 
on HLA hull in hover. T /S R * 395 N/m* (8.26 psf). 







Solid vectors are in-ground effect 
Dashed vectors are out-of-ground effect 
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Figure 18.- Rotor- induced velocities on hull at 
a ■ 5 m (16 ft) for HLA in axial flow at 
15.4 m/s (30 knots) in- and out- 
of-ground effect. 
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Figure 21.- Rotor-induced velocities on HLA hull 
in- and out-of-ground effect, a - 0 # p * -90°. 
x - -15m (-48 ft). 
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